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Abstract 


By analyzing 2.93 #7! of data taken at the (3770) resonance peak with the BESIII detector, we measure the 
branching fractions for the hadronic decays Dt > K2K2Kt, Dt > ККОт+, D? > КОК and D? > 
КОКОКО. They are determined to be B(D* > K3K3K*) = (2.54 + 0.05stat. + 0.12sys.) x 107%, B(Dt > 
КОКО) = (2.70 + 0.0544. + 0.12sys.) x 1073, BUD? — KORE) = (1.67 + 0.11stat. + 0.115) x 1074 and 
B(D? + КОКК) = (7.21 £0.33stat. + 0.44554.) x 1074, where the second one is measured for the first time and 
the others are measured with significantly improved precision over the previous measurements. 


Keywords: BESIII, D? and D* mesons, Hadronic decays, Branching fractions. 


1. Introduction 


Hadronic decays of D mesons open a window to 
probe for the physics mechanisms in charmed meson 
decays, e.g., CP violation, D? D? mixing and SU(3) 
symmetry breaking effects. Since the discovery of D 
mesons in 1976, the hadronic decays of D mesons have 
been extensively investigated [1]. However, the existing 
measurements of the D hadronic decays containing at 
least two Kg mesons in the final state are still very poor 
due to limited statistics [1]. 

In this Letter, we report the measurements of the 
branching fractions for the hadronic decays Dt — 
K9K8-*, D? + КОКО, Dt > КОКОК+ and 
D? — КОКК. Throughout this Letter, charged 
conjugate modes are implied. These decays have sim- 
pler event topologies and suffer less from combinato- 
rial backgrounds than other decay modes containing 
two K 2 in the final state. The comprehensive or im- 
proved measurements of three-body decays will benefit 
the understanding of the interplay between the weak and 
strong interactions in multibody decays where theoret- 
ical predictions are poorer than two-body decays. The 
improved measurements of two-body decays can serve 
to better explore the contributions of W-exchange dia- 
grams and final-state interactions [2, 3, 4, 5], as well as 
SU(3)-flavor symmetry breaking effects [6, 7, 8, 9, 10] 
in D meson decays. In addition, these measurements 
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will also help to improve background estimations in the 
precision measurements of D and B meson decays. 

The data sample used for this analysis, which has an 
integrated luminosity of 2.93 fb ^! [11], was taken at the 
~(3770) resonance peak with the BESIII detector [12]. 
The D? D? and D+ D~ pairs produced in (3770) de- 
cay provide cleaner D? and D* meson samples than 
those used in previous studies at ARGUS [13, 14], 
CLEO [15, 16] and FOCUS [17]. To optimize the 
precision for these measurements, we use a single-tag 
method, in which either a D or D is reconstructed in an 
event. We combine the yields measured with previously 
reported values of the cross sections for e* e- — D? D? 
and D* D- at the (3770) resonance peak [18]. 


2. BESIII detector and Monte Carlo simulation 


The BESIII detector is a magnetic spectrometer that 
operates at the BEPCII collider. It has a cylindrical ge- 
ometry with a solid-angle coverage of 93% of 4m. It 
consists of several main components. A 43-layer main 
drift chamber (MDC) surrounding the beam pipe per- 
forms precise determinations of charged particle trajec- 
tories and measures the specific ionization (dE /dx) for 
charged particle identification (PID). An array of time- 
of-flight counters (TOF) is located outside the MDC and 
provides additional PID information. A СЅКТІ) elec- 
tromagnetic calorimeter (EMC) surrounds the TOF and 
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is used to measure the energies of photons and elec- 
trons. A solenoidal superconducting magnet outside 
the EMC provides a 1 T magnetic field in the central 
tracking region of the detector. The iron flux return of 
the magnet is instrumented with 1272 m? of resistive 
plate muon counters (MUC) arranged in nine layers in 
the barrel and eight layers in the endcaps for identifica- 
tion of muons with momentum greater than 0.5 GeV/c. 
More details about the BESIII detector are described in 
Ref. [12]. 

A GEANTA-based [19] Monte Carlo (MC) simu- 
lation software package, which includes the geometric 
description and response of the detector, is used to de- 
termine the detection efficiency and to estimate back- 
ground for each decay mode. An inclusive MC sam- 
ple, which includes the D°D°, D+ D- and non-DD 
decays of the 1» (3770), initial-state-radiation (ISR) pro- 
duction of the (3686) and J/w, the ete” — qq 
(q = u, d, s) continuum process, the Bhabha scatter- 
ing events, the di-muon events and the di-tau events, is 
produced at \/s = 3.773 GeV. The equivalent luminos- 
ity of the MC sample is ten times of data. The (3770) 
decays are generated by the MC generator KKMC [20], 
which incorporates both ISR effects [21] and final-state- 
radiation (FSR) effects [22]. Known decay modes are 
generated using EvtGen [23] with input branching frac- 
tions from the Particle Data Group (PDG) [1]. Unmea- 
sured decays are generated using LundCharm [24]. 


3. Data analysis 


All charged tracks used in this analysis are required 
to be within a polar-angle (0) range of |cos 0| < 0.93. 
The good charged tracks, except when used to recon- 
struct КО mesons, are required to originate within an 
interaction region defined by Vzy < 1.0 cm and У, < 
10.0 cm, where Vzy and V; are the distances of clos- 
est approach of the reconstructed track to the interaction 
point (IP) perpendicular to (xy) and along (2) the beam 
direction. 

The charged kaons and pions are identified by the 
dE/dx and TOF measurements. The combined confi- 
dence levels for pion and kaon hypotheses (СГ. and 
C Lx) are calculated, respectively. The charged track 
is identified as kaon (pion) if Ск > CL, (CLr > 
CL x) is satisfied. 

K candidate mesons are reconstructed through the 
mm decay mode. Charged pions used in Ke can- 
didates mesons are required to satisfy V, < 20.0 cm. 
The two oppositely charged tracks are assumed to be a 
тл pair without PID requirements. To reconstruct 
K 2, the тїт combination is constrained to have a 


common vertex. The candidate is accepted if it has an 
invariant mass M,+,- within 12 MeV/c? of the Kg 
nominal mass [1] and satisfies L/or, > 2, where L is 
the measured flight distance and су, is its uncertainty. 

To identify D candidates, we use two selec- 
tion variables, the energy difference AE = Ep — 
Ebpeam and the beam-energy-constrained mass Mpc = 
МЕ? м/с — |Po|?/c?, where Ebeam is the beam en- 
ergy and Ep and pp are the energy and momentum 
of the D candidate in the e+ e^ center-of-mass system. 
For each signal decay mode, only the combination with 
the minimum |AF| is kept in events where more than 
one candidate passes the selection requirements. Mode- 
dependent A E cuts are determined separately for data 
and MC based on fits to the respective AE distributions. 
These are set at +30, where c is the A E resolution (Ta- 
ble 1). 


The combinatorial т * 7 | ко pairs with invari- 


ant mass in КО signal region may also satisfy the K2 
selection criteria and contribute peaking background 
around the D mass in the Мвс distribution. This peak- 
ing background is estimated with events in the К 0 side- 
band region, defined as 0.020 < |M,+,- — Мко| < 
0.044 GeV/c?. Figure 1(a) shows the comparison of 
the M,,+,- distribution for D? +> КОКО candidates 
in data with the corresponding distribution for the in- 
clusive MC. In the figure, the solid (dashed) arrows de- 
lineate the К р signal (sideband) regions. 


In the analyses of the D? + КОКО, Dt — 
K9KQK* and КОКОлт decays, two-dimensional 
(2D) signal and sideband regions are defined. Fig- 
ure 1(b) shows the distribution of M,+,-(1) versus 
M,+x-(2) for the D? — КОКО candidate events in 
data. The solid box, in which both of the 7+ 27 com- 
binations lie in the КО signal regions, denotes the 2D 
signal region. The dot-dashed (dashed) boxes indicate 
the 2D sideband 1 (2) regions, in which one (two) of 
the л+л combinations lie in the K E sideband regions 
and the others аге in the K$ signal region. For the 
D? + КЎК КО decay M,«,-(1 versus Mz, (2) 
versus Мух+х-(зу of the candidate events in data is 
shown in Fig. 1 (c). The region in which all three 7* x^ 
combinations lie in the К E signal regions is taken as the 
three-dimensional (3D) signal region. The 3D sideband 
i (i = 1,2,3) regions denote those in which i of the 
three mtr” pairs lie in ће K$ sideband regions and 
the rest are located in the К 3 signal regions. 

The resulting Mpc distributions of the accepted 
candidate events in the 2D or 3D signal region, sideband 
1 region and sideband 2 region are shown in the sub- 
figures of the first, second and third rows of Fig. 2, re- 


TABLE 1: AF requirements (in MeV) for data and MC samples. 


Decay modes Data MC 
D? + КУКУК“ (—17,419) (—16,+16) 
Dt — K3K8n*  (—-175,447) (17, +16) 
D? + кк (—19, +17) (—17, +14) 
D? КОКОКО  (—14,416) (—13, +13) 
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FIG. 1: (a) Comparison of the M +p- distributions of the Do > K оК Е candidate events between data (dots with error bars) and inclusive MC 
(histogram). The pairs of the solid (dashed) arrows denote the К Е signal (sideband) regions. (b) Distribution of М. +. (1) Versus М.а (2) 


for ће D? — K9 Ke candidate events in data. (c) Distribution of M,.+,,—(1) versus M,,+,,—(2) versus M4. — (зу for the р? + кококо 
candidate events in data. In these figures, all selection criteria have been imposed except for the К n mass requirement and Мвс is required to be 
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within 5 MeV/c? around the D nominal mass [1]. 


FIG. 2: Fits to the Myc distributions of the (a) Dt  KOK9K*,(b) Dt  КОКОт+, (с) D? > КОКО and (d) D9 > КОКОКО 
candidate events. The dots with error bars are data, the solid curves are the total fits, and the dashed curves are the fitted backgrounds. The 
first, second and third rows correspond to the fits to the candidate events in the 2D or 3D signal region, sideband 1 region and sideband 2 region, 
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spectively. By fitting these Мвс distributions as shown 
in Fig. 2, we obtain the fitted yields of D signal in the 
2D or 3D signal region, sideband 1 region and sideband 
2 region, IN ko sig: Nsb1, Nsb2, which are given in Ta- 
ble 2. In the fits, the D signal is modeled by a MC- 
simulated shape convoluted with a Gaussian function 
with free parameters accounting for the difference of 
detector resolution between data and MC. The combi- 
natorial backgrounds are described by an ARGUS func- 
tion [25] with an endpoint of 1.8865 GeV/c”. In the 
Мвс fits for the 2D or 3D sideband events, the param- 
eters of the convoluted Gaussian function are fixed at 
the values determined for the signal region. For the 
D? — КОКОКО decays, the peaking backgrounds 
from sideband 3 region are negligible since few events 
survive. 


In this analysis, the combinatorial background in the 
M,,+,- distribution are assumed to be flat, which im- 
plies that the ratio of background yields between the 
K? signal and sideband regions is 0.5. Thus, the net 
numbers of the D? + КОКО, Dt > КОКК and 
K&K&n* decays can be calculated by 


1 1 
Мое = Мко Мзь1 + 12 = N? өг (1) 


sig 9 


and the net number of the D? + K} КО КО decays сап 
be calculated by 


1 ‚1 1 b 
Мое = Мкозш shi T т\зъ? 8 Мыз Nother 
(2) 
where Акош and Nsp; are D signal yields from the fit 


in the 2D or 3D signal regions and sideband i regions, 
respectively. Л н is the normalized number of resid- 
ual peaking background. For ће Dt > КОКОК+, 
D* + K2K8n* and D? —^ КОКОКО decays, the 
residual peaking background is mainly from the events 
of Dt 2 KRKP Kt, Dt > KOK9n* and D? > 
КОКК? versus D (D?) + КОХ (X = any pos- 
sible particle combination). This kind of background 
peaks around the nominal D mass [1] when the Ke 
from а D~(D°) decay has momentum similar to that 
of a K? produced in D* (D?) decay. These peaking 
backgrounds cannot be modeled by the events from the 
2D or 3D sideband region and are estimated by analyz- 
ing the inclusive MC sample. The measured values of 
NP, .. and Naet are given in Table 2. 


other 


4. Branching fractions 


The branching fraction for the hadronic decay 
D+) — f is determined by 


Nnet 


B(Dt® — one 0 00. 
( +f) 2-0p«p- (popo): L: e’ 


(3) 


where Ме is the net number of DH a f decays in 
data, є is the detection efficiency including the branch- 
ing fraction of КО — m7, С is the integrated lu- 
minosity of data [11] and ср+р- (po po) is the D D^ 
(D? D?) cross section at the (3770) resonance peak. 

The detection efficiencies are determined by an- 
alyzing the inclusive MC sample. In this sample, 
the signal MC events for Dt — K&KQn* are 
produced as a mixed sample containing 9096 of the 
D+ — K2K*(892)*, K*(892)* — Km" decays 
and 10% of the direct three-body decay in phase space 
D* — K$K$m'. The signal MC events for D+ + 
Кокок, D? — КОКО and КОКОКО are pro- 
duced using a phase-space model. Detailed studies 
show that the momentum and polar-angle distributions 
of the daughter particles in data are well modeled by 
the MC simulation for each decay mode. By analyzing 
the inclusive MC sample with the same analysis pro- 
cedure applied to the data (including the Мвс fits and 
the calculation of the net signal yields), we obtain the 
net number of D mesons observed for each decay. The 
detection efficiency e is obtained by dividing the net D 
signal by the total number of signal events, taking into 
account the efficiency correction discussed in Sect. 5. 

Inserting the numbers of Nnet, є, £, as well as 
0p-p- = (2.882 + 0.018stat. + 0.042,,,) nb or 
opopo = (3.607 + 0.017stat. + 0.056,,, ) nb quoted 
from Ref. [18] into Eq. (3), we obtain the branching 
fraction for each decay, as listed in Table 2, where the 
uncertainties are statistical only. 


5. Systematic uncertainty 


Table 3 shows the systematic uncertainties in the 
branching fraction measurements. Each of them, es- 
timated relative to the measured branching fraction, is 
discussed below. 


* MC statistics: The uncertainties due to the limited 
MC statistics are 0.5%, 0.4%, 1.8% and 1.3% for 
Dt — KLK} K+, Dt > КоКОт+, D9 + 
K} K} and D? + КОКОКО, respectively. 


* Luminosity of data: The uncertainty in the quoted 
integrated luminosity of data is 0.5% [11]. 


TABLE 2: Input quantities and results for the determination of the branching fractions as described in the text. The uncertainties are statistical only. 


Decay modes Ngo sig Nap Ngp2 
D? + КУКУК" 3616+ 66 97+19 6+ 8 
D* > K8K8nt 5643 +88 1464+68 69+19 
D? коко 888 +36 626431 3+ 6 
D? КОКОКО 622427 24+ 8 14+ 6 


DD cross section: The uncertainties of the 
quoted D* D- and 00 00 cross sections are 
1.6% [18]. 


В(КО — ntr): The uncertainty of the quoted 
branching fraction for K? — път is 0.196 [1]. 


K reconstruction: The K reconstruction ef- 
ficiency has been studied as a function of mo- 
mentum by using the control samples //ф — 
K*(892)* K* and J/$ > ФК Кт. Small 
data-MC efficiency differences are found and pre- 
sented in Ref. [26]. To correct the Ke recon- 
struction efficiency, a piecewise fit to these dif- 
ferences as a function of K} momentum is per- 
formed. For the efficiencies of detecting the de- 
cays Dt > КОКК, Dt > KOKOm*, 
D? — K$K$ and D? — КОКОКО, the mo- 
mentum weighted differences associated with К п 
reconstruction between data and MC are deter- 
mined to be (+3.9 + 1.9)%, (4-3.0 + 1.4)%, 
(+1.8 + 0.8)% and (+5.9 + 2.8)%, respectively, 
where the uncertainties are statistical. These cor- 
rections are applied to the detection efficiencies, 
after which only the statistical uncertainties of the 
differences are retained. On average, the residual 
uncertainty for each КО is no more than 1.0%. 
Furthermore, the difference of the momentum- 
weighted efficiencies between data and MC from 
the different fits, which is 1.0% per K 0. is in- 
cluded as an additional uncertainty. Finally, we 
assign 1.5% per K$ as the systematic uncertainty 
for the reconstruction efficiency . 


Tracking [PID] for K * (1^): The tracking [PID] 
efficiencies for K* and т* are investigated us- 
ing doubly tagged DD hadronic events. The dif- 
ference of momentum weighted efficiencies be- 
tween data and MC of the tracking [PID] are de- 
termined to be (+2.1 + 0.4)% [(—0.3 + 0.1)%] 
for ће K* in ће Dt + КОКОК+ decay and 
(+0.4 + 0.3)% [(—0.3 + 0.1)%] for the т in 
the Dt — КО Котт decay, where the uncertain- 
ties are statistical. After correcting the detection 
efficiencies by these differences, we take 0.5% 


Nsb3 


Nother Nue e (96) B (x105 
1822 3551467 8271x004 254405 
3143 4897494 10.7240.04 27.040.5 
0 576 + 39 16.28 +0.30 1.67 0.11 
16+3 597427 3.92+40.05 17.21 0.33 


[0.5%] as the systematic uncertainties in tracking 
[PID] for ће K+ and тт, respectively. 


Мвс fit: In order to estimate the systematic un- 
certainty associated with the Мвс fit, we re- 
peat the measurements by varying the fit range 
((1.8415, 1.8865) GeV/c?), signal shape (with 
different MC matching requirements) and end- 
point of the ARGUS function (+0.2 MeV/c?). 
Quadratically summing the changes of the 
branching fractions yields 2.196, 1.096, 4.296 and 
2.7% for Dt + K8K9K^, D* + ККОт+, 
D? + КК and D? + КОККО, which are 
assigned as the relevant systematic uncertainties. 


AE requirement: 'To investigate the systematic 
uncertainty associated with the A E requirement, 
we repeat the measurements using alternative AE 
requirements of +(4,5,6) times the resolution 
around the AE peaks. The maximum changes 
of the branching fractions, 2.0%, 1.5%, 2.0% and 
1.5% for D* + K2K2K*, D* + КОКОт+, 
D? + КОКО and D? + КОКОКО, are taken 
as the associated systematic uncertainties. 


Normalization of peaking backgrounds: In the 
nominal analysis, the normalization factor for the 
peaking backgrounds, which is the ratio of back- 
ground yields between the K? signal and side- 
band regions, has been assumed to be 0.5. The 
branching fractions are recalculated with alter- 
native normalization factors determined by MC 
simulation. The corresponding changes on the 
branching fractions, 0.5%, 1.4%, 2.4% and 0.7% 
for Dt > KÌKĮK+*, Dt > KQK$m-', 
D? — КК and D? > КОКОКО, are as- 
signed as the systematic uncertainties associated 
with the peaking background (PBKG) normal- 
ization. On the other hand, the uncertainties 
of the residual peaking backgrounds are domi- 
nated by the uncertainties of the input branch- 
ing fractions for D~(D°) — КОХ, which con- 
tribute additional uncertainties of 0.1%, 0.1% and 
0.4% for the measured branching fractions for 


Dt + KÜKÜK*, Dt > К Кол+ and D? > 
KÌ K} K9, respectively. 


K 9 sideband: To evaluate the systematic uncer- 
tainty due to the choice of K E sideband region, 
we remeasure the branching fractions after shift- 
ing the КО sideband by +2 MeV/c?. The cor- 
responding maximum changes in the branching 
fraction, which are 0.596, 0.596, 2.096 and 1.096 
for Dt + KGKSK*,D* > КОКОт+, Dh > 
K9K$ and D? — КОКОКО, respectively, are 
taken as the systematic uncertainties. 


* MC modeling: For the three-body decays, we ex- 
amine the reweighted detection efficiencies by in- 
cluding the possible sub-resonances ао (980) and 
fo(980) in the signal MC samples. The maxi- 
mum change of the reweighted detection efficien- 
cies, 1.0%, is taken as the systematic uncertainty 
in MC modeling. 


Adding all of above systematic uncertainties in 
quadrature, we obtain the total systematic uncertainties 
of 4.7%, 4.4%, 6.8% and 6.1% for Dt + КОКК, 
D+ > ККт+, D D + КК and D > 
KÌ KÌ K9, respectively. 


6. Summary 


In summary, by analyzing 2.93 fb^! of data col- 
lected at ys = 3.773 GeV with the BESIII de- 
tector, we measure the branching fractions for the 
hadronic decays Dt + K8K9K*, Dt > КО КОт+, 
D? — ККЗ and D? > КОКОКО using a single- 
tag method. Table 4 presents the comparisons of 
the measured branching fractions with the PDG val- 
ues [1]. The branching fraction for Dt > KS K&z* 
is measured for the first time and the others are con- 
sistent with previous measurements, but with much 
improved precision. We also determine the branch- 
ing fraction ratios B(D* —^ K&K&K*)/B(Dt > 
KeKent) = 0.941 + 0.025... + 0.040sys. and 
B(D? ^ K3K9)/B(D® > K8K8K$) = 0.232 + 
0.019stat. Е 0.0165,,., in which the systematic uncer- 
tainties in the D+ D~ (or D? D^) cross section, the in- 
tegrated luminosity of data, as well as the reconstruction 
efficiencies and the branching fractions of the two Kg 
mesons cancel. The results in this analysis provide help- 
ful experimental data to probe for the interplay between 
the weak and strong interactions in charmed meson de- 
cay [2, 3, 4, 5]. In addition, the measured branching 
fraction for the two-body decay D? —^ КОКО can also 


help to understand SU(3)-flavor symmetry breaking ef- 
fects in D meson decays [6, 7, 8, 9, 10]. 
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